#!/usr/bin/env python3
# -*- coding: utf-8 -*-
"""
Created on Mon Aug 23 22:14:19 2021

@author: liqingsimac
"""

## Example 7-5-1
import numpy as np
import matplotlib.pyplot as plt
import scipy.special as spl 

x = np.linspace(0,20,101)
y0 = spl.jv(0,x)
y1 = spl.jv(1,x)
y2 = spl.jv(2,x)

fig=plt.figure()
ax=fig.add_subplot(111)
ax.plot(x, y0, 'r-', label='n=0')
ax.plot(x, y1, 'g-', label='n=1')
ax.plot(x, y2, 'b-', label='n=2')
ax.legend(loc='best')
fig.savefig('ode-example-7-5-1.png')

